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In recent letter [Phys. Rev. Lett 105, 150606 (2010)], the surface-mediated diffusion problem 
is theoretically discussed, and interesting results have been obtained. However, for more general 
cases, the ansatz of solutions of the diffusion equation, which is the starting point of their analysis, 
might not be appropriate. In this comment, suggested ansatz and corresponding methods will be 
presented. 
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For the surface-mediated process studied in in 
which a molecule diffuses in a spherical confining do- 
main of radius R and alternates phases of boundary dif- 
fusion (with diffusion coefficient Di) and phases of bulk 
diffusion (with diffusion coefficient D2), the mean first- 
passage (MFPT) at the target, i.e., parts of the boundary 
r = i?: e [0, e] U [27r — e, 2-k] , satisfies the following equa- 
tion 


t'l{e) - aXtiiO) = -(7[l + \t2{R-a,e)], 


(1) 


for e e [e,27r - e]. Where ti{e) is the MFPT starting 
from the boundary defined by angle 6, A is dissociation 
rate from the boundary to bulk, a is the distance that the 
molecule is assumed to be ejected once it is dissociated 
from the boundary, and a = R^/Di. t2{r, 9) is the MFPT 
staring from point (r, 0) , which satisfies 


D2{d^ + dr/r + d^/r^)t2{r, 9) = -1. 


(2) 


The boundary conditions of Eqs. ([T]) ([2]) are t2{R,0) = 
i 1 (6') and i 1 (6*) = for e [0, e] U [27r - e, 27r] respectively. 

In [H and also in their more recent work [3], the fol- 
lowing ansatz is assumed to be the solutions of Eq. ^ 


t2{r,0) = ao 


00 


— -f >^a„r"cos(n0), (3) 


n=l 


where a„ are constants to be determined. Based on this 
ansatz, interesting results have been drawn in [1]. How- 
ever, mathematically, it can be easily shown that a more 
general form of the solutions of Eq. ([2]) is 


most of the physical or biophysical processes, it is rea- 
sonable to assume that t2{0,0) < 00. Therefore, the 
general form 1^ can be simplified to 


t2{r..0) =ao(l + 7^) 


r"[an cos{n9) + /3n sm{n6)], 


(5) 


n=l 


In a general aspect, the ansatz ([SJ might be more appro- 
priate than ansata ([3]). 

To obtain the MFPT (ti), in Eq. (P) is solved 
approximately by perturbative expansion, and then an- 
other approximate method is developed in Here, I 
will suggest another method. For the general form ([5]) 
of t2{r,d), the solution ^1(6*) of Eq. ([1]) can be obtained 
explicitly (we assume R and A ^ 0), 


11(0) — C\e^ -f C2e 


-f ao(l +761) - 


{R-a) 


1 
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(6) 


(i? — a)" 

„ , . J an cos(n6') + /3„ sin(n6')], 
^ — ' n^/crA -f- 1 

n— 1 ' 


where constants ci , C2 are determined by boundary con- 
ditions ti(e) ~ ii(27r — e) = 0, and aQ,an,l3n can be 
obtained by t2{R,0) = ti{e) and Eqs. (P and di]). 

If a = i? or A = 0, then Eq. fl]) is reduced to t'{{0) = 
— (7. Its solution is 


ti{d) = a{2TT - e-e){e - e)/2. 


(7) 


t2{r,e) =ao(l + A*lnr)(l+7(?) 
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(4) 


r'^[an cos{n6) + (3n sin(n6')]. 


where Q;„,/3„,/i,7 are constants to be determined. For 


Then the MFPT (ti) can be easily obtained by integral 
"'tide. 

It should be pointed out that the above discussion is 
only valid for A is constant. In real physical process, A 
will decrease with a, and A cxd when a ^ 0. Therefore 
the formulation ^ is only an instructive approximation. 
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